Astronomy & Astrophysics manuscript no. 9897 
July 28, 2008 



© ESO 2008 



Multiplicity of young stars in and around R Corona Australis' 

R. Kohler^'^, R. Neuhauser^ S. Kiamer^ Ch. Leinert^, T. Ott"^, and A. Eckart^ 



00 

o 
o 

(N 



00 

6 



> 

cn 
a^ 
m 
^. 

o 

oo 
O 



X 



ZAH Landessternwarte, Konigstuhl, 69117 Heidelberg, Germany, e-mail: r.koehlerOl sw.uni-heidelberg.de 

Max-Planck-Institut fiir Astronomie, Konigstuhl 17, 691 17 Heidelberg, Germany 

Astrophysikalisches Institut und Universitats-Sternwarte, Friedrich-Schiller-Universitat Jena, Schillergasschen 2-3, 07745 Jena, 

Germany 

MPI fiir extraterrestrische Physik, Postfach 1312, D-85741 Garching, Germany 

I. Physikalisches Institut, Universitat zu Koln, Zulpicher StraBe 77, 50937 Koln, Germany 



Received April 3., 2008; accepted June 26., 2008 



ABSTRACT 



Context. In star-forming regions like Taurus-Auriga, it has been found that most young stars are born as multiples, which theories 

for star formation should definitely take into account. The R CrA star-forming region has a small dark cloud with quite a number of 

protostars, T Tauri stars, and some Herbig Ae/Be stars, plus a number of weak-line T Tauri stars around the cloud found by ROSAT 

follow-up observations. 

Aims. We would like to detect multiples among the young stars in and around the R CrA cloud in order to investigate multiplicity in 

this region. 

Methods. We performed interferometric and imaging observations with the speckle camera SHARP I at the ESO 3.5 m NTT and 

adaptive optics observation with ADONIS at the ESO 3.6 m telescope, all in the near-infrared bands JHK obtained in the years 1995, 

2000, and 2001. 

Results. We found 13 new binaries among the young stars in CrA between 0.13 arcsec (the diffraction limit) and 6 arcsec (set as an 

upper separation limit to avoid contamination by chance alignments). While most multiples in CrA are binaries, there are also one 

quadruple (TY CrA), and one triple (HR 7170) which may form a quintuple together with the binary HR 7169. One of the newly 

detected companions with a large magnitude difference found near the M3-5 type T Tauri star [MR 81] Ha 17 could be a brown dwarf 

or an infrared companion with an edge-on disk. Among seven Herbig Ae/Be stars in CrA, six are multiple. 

Conclusions. The multiplicity frequency in CrA is as high as in similar star forming regions. By comparing with the period distribution 

of main-sequence stars and extrapolating to separations not probed in this survey, we conclude that the companion-star frequency is 

(95 ± 23) %; i.e. the average number of companions per primary is 0.95. 

Key words. Stars: pre-main-sequence - Binaries: close - Infrared: stars - Instrumentation: high angular resolution - Surveys 



1. Introduction 

The star-forming region Corona Australis (southern crown), ab- 
breviated CrA or R CrA, is today known as one of the nearest re- 
, gions of ongoing and/or recent intermediate- and low-mass star 
■ formation. The dark cloud around R CrA has an extinction of up 
to Ay ~ 45 mag. The age ranges between < 1 Myrs for the pro - 
tost ars to 9 ± 4 Myrs for some T Tauri stars (iJames et al.ll2006h . 
See Neuhauser & ForbrichI ( |2008|) for a recent review. 

ICasev et all (Il998l) determined the distance towards the 
eclipsing double-lined spectroscopic binary (SB) TY CrA 
to be 129 + 11 pc from their orbit solution. Because this 
is consistent with other ea rlier estimates (see discussion in 
iNeuhauser & Forbrichl2008h and, at the same time, the best cur- 
rent estimate, we use ~ 130 pc as the distance for the young stars 

in CrA. 

INeuhauser & ForbrichI (|2008|) give a list of all known opti- 
cally visible members compiled from various sources. Only a 
few multiples have been published among the CrA members so 
far; see Table 1 for a listing. Hence, while six out of seven Herbig 
Ae/Be stars are known or suspected to be multiple, only five 
T Tauri stars (TTS) were known to be multiple, a low number 



* Based on observations obtained at the European Southern 
Observatory, La Silla, proposal numbers 55.E-0968, 65.1-0086, and 
67.C-0213 



compared to other star- forming regions. In a ddition, the brow n 
dwarf member Denis 1859 is a known binary dBouv et al.ll2004l) . 
Therefore, we performed this new homogeneous multiplic- 
ity survey of most CrA members with the infrared speckle cam- 
era SHARP I and the infrared adaptive optics (AO) instrument 
ADONIS. We explain our observations and the data reduction in 
Sect. 2, list all results in Sect. 3, compare our observations of 
multiples with previous observations in Sect. 4, and conclude in 
Sect. 5. 

2. Observations and data reduction 

We have observed the sample of 49 optically visible young mem- 
bers of the R CrA association listed in Tabled This list i s com - 
pile d fromlGlass & Pens tonI (1 1 975l). lMarraco & RvdgrenI (1 1 98 lb. 
Her big & Belli d 19881). I W alter et all d 19971) . as well as lNeuhauseii 
(1997) and'Neuhau.sereLalj (|200§. 

The majority of our targets (Table |2]i were observed with 
the speckle interferometry method during two observing runs in 
July 1995 and July 2001 at the European Southern Observatory 
(ESO) 3.5 m New Technology Telescope (NTT) on La Silla, 
Chile. We used the SHARP I camera (System for High 
Angular Resolution Pi ctures) of the Max-P lanck-Institut for 
Extraterrestrial Physics dHofmann et alj[l992h . All observations 
were done in the K-band at 2.2 /im; some of the stars, where a 
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Table 1. Previously known multiples among optically visible young stars in and around CrA (c,d) 



Designation 



Separation 
[arc sec] 



Position 
Angle [°] 



Brightness ratio (, 
H 



a) 



K 



Epoch 
day/month/year 



Reference 



SCrA 



1 


135 


1 


134 


1.41 ±0.06 


157 ±2 


1.37 ±0.02 


147 ±2 


1.3 


149 


.346 ± 0.009 


144.5 ± 0.2 


1.34 ±0.04 


147 



0.45 
0.49 ± 0.03 



0.41 
0.49 ± 0.04 





04/07/1942 


(1) 
(15) 


0.3 ± 0.2 (b) 


12/04/1995 


(2) 




1981.25 


(3) 




April 1991 


(4) 




25/08/1996 


(5) 


0.39 


06-08/06/1987 


(6) 


0.52 ± 0.05 


01/07/1996 


(7) 



TY CrA B-C 

A-BC 

A-D 


eclipsing spectroscopic binary 
4.1 138 
0.294 ±0.007 188.5 ±1 31/03/2002 


(13,14) 
(15) 
(8) 


RCrA 


spectro-astrometric binary according to Takami et al. 2003, but no binarity detected in Bailey 1998 




TCrA 


> 0.076 ±0.005 273.0 ±1.4 26/08/1996 

> 0.140 ±0.009 277.7 ±1.3 26/06/1997 


(5) 
(5) 



CrA 


2.1 ±0.1 


48 ±2 




2.08 ± 0.025 


46.5 ±0.12 




1.9 


44 




2.10 


44 



0.98 
0.1 ±0.02 



0.26 
0.3 ± 0.03 



0.96 ± 0.02 (b) 



0.12 
0.95 ± 0.06 



12/04/1995 


(2) 


1993/1994 


(9) 


April 1991 


(4) 


06-08/06/1987 


(6) 


01/07/1996 


(7) 



CrAPMS 3 


4.5 
4.66 


58 
58 


0.12 
0.14 ±0.01 


0.13 0.15 
0.19 
0.15 ±0.01 0.17 ±0.02 


April 1991 
06-08/06/1987 

01/07/1996 


(4) 
(6) 
(10) 

(7) 


CrAPMS 6 


3.71 


207.26 


0.95 


0.95 0.89 


1997 


(10) 


RXJ1857.5-3732 


<4 










(11) 


DENIS-PJ 

185950-370632.9 


0.065 ± 0.001 
0.057 ± 0.0005 
0.066 ± 0.003 
0.059 ± 0.003 


279.2 ±0.1 
279.1 ±0.1 
283.8 ± 1.2 
271.8 ±1.2 






24/09/2002 
24/09/2002 
12/09/2000 
12/09/2000 


(12) 
(12) 
(12) 
(12) 


HR7169 


spectroscopic binary (WDS. Worlev & Doualass 


1997). 13 arc sec CPM pair with HR 7170 




HR7170 


spectroscopic binary (WDS, Worlev & Douglass 1997), 13 arc sec CPM pair with HR 7169 




HD 176386 


^4 


~ 137 








(16) 



Remarks: (a) Always < 1 here by definition, (b) For consistency, we give here the reciprocal value of the flux ratio given by Ghez et al. (c) 
RXJ 1846.7-3636 was listed as ~ 8" binary in (11), i.e. with separation above t he 6" upper limit used h ere, hence regarded as two separate stars 

here, (d) The very low-mass star or brown dwarf CrA 444 may also be a binary (Lopez Marti et al. 2005). 

References: (1) Jov & van Biesbroeck ( 1941), (2)jGhez et al. ( 1997), (3) Baier et al. ( 1985), (4) Reipurth & Zinneckej (Il993h . ( 5)lBailevllll99^ 
(6)[Chelli et al. (1995), (7) Prato et al. (2003), (8) Chauvin et al. (2003), (9) Ageoraes et al. (1997), (10) Walter et al. ( 1997), (11) Neuhauser et al] 
( 2000) , (12) .Bouvetal..c2004.) . (13) .Gaposchkin & Greenstein H936.) . (14) .Casev et al..H998.) , tl5) .Proust et al..tl981.) (16) .Dommanget & Nv3 
( 20021) 



companion candidate was detected, were also observed in the H- 
and J-bands at 1 .2 fim and 1 .6 fim, respectively. For most targets, 
integration times of 0.5 s per frame were used. For the bright 
stars R and T CrA, the time was reduced to 0.2 s to avoid sat- 
urating the detector. On each target, 500 frames were taken in 
1995, or 600 frames in 2001, giving a total integration time of 
100 s for R and T CrA, and 250 - 300 s for the other targets. 
To allow background subtraction and bad pixel correction, the 
telescope was moved after half of the frames had been taken, to 
position the target at a different spot on the detector 

A few targets were observed in June 2000 with the AO sys- 
tem ADONIS and the SHARP II camera at the ESO 3.6 m tele- 
scope on La Silla. We used the K-band filter of this instrument, 
which has a central wavelength of 2.177/im. The observing 
strategy was similar to that used for the speckle observations, 
i.e. we took many frames with short integration times. This al- 
lows us to use the same programs and algorithms for data reduc- 
tion. Since the AO system corrects the atmospheric turbulence, 
somewhat longer integration times of 1 - 2 s per frame could be 
used. The targets were observed at 4 different positions on the 



detector. In total, 120 - 240 frames were recorded to obtain a to- 
tal integration time of 240 s for each target. Both SHARP I and 
ADONIS/SHARP II were used in a configuration with a field- 
of- view of about 12" x 12". 

The detectors of the SHARP cameras use separate read-out 
electronics for the four quadrants, which leads to discontinuities 
at the quadrant borders. To avoid distortions in the images of our 
stars, we positioned the targets in the centers of one quadrant. 
This implies that companions at separations larger than about 3" 
might be outside the field-of-view. To make sure we did not miss 
any co mpanions, we s earched the 2MASS point source catalog 
tSkrutskie et al.ll2006l) for sources near our targets. We did not 
find any additional companions within 6". 

Although speckle i nterferometry can be considered by now 
a standard technique ( Leinert|[l994l) . no program for speckle 
data reduction was publicly available at the time this survey was 
started. Therefore we used our speckle program, which was al- 
ready used for binary survey s in a number of other star-forming 
regions, e. g. Taurus-Auriga ( Kohler & Lein ert 1998 ), Scorpius - 
Centaui-us dKohler et al.ll2000h . and Chamaeleon (,K5hleiil2001h . 
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Table 2. Stars observed in this work 



No. Designation 



Q'2000 



Date(s) of Observation(s) lustrum. Filter(s) Type (a) SpTy Note 



1 


SCrA 


19:01:08.6 


-36:57:20 


12 July 1995, 06 July 2001 


SHARP 


JHK 


c 


K6 


2 


TYCrA 


19:01:40.8 


-36:52:34 


12 July 1995 


SHARP 


K 


H 


B9 


3 


RCrA 


19:01:53.6 


-36:57:08 


12 July 1995 


SHARP 


K 


H 


A5 


4 


DGCrA 


19:01:55.2 


-37:23:41 


12 July 1995 


SHARP 


K 


c 


KO 


5 


TCrA 


19:01:58.8 


-36:57:50 


12 July 1995 


SHARP 


K 


H 


FO 


6 


VVCrA 


19:03:06.7 


-37:12:50 


12 July 1995, 06 July 2001 


SHARP 


JHK 


c 


Kl 


7 


[MR81] Ha 10 


18:58:51.8 


-37:19:23 


12 July 1995 


SHARP 


K 


n 


K (b) 


8 


[MR81] Ho- 6 


19:00:01.6 


-36:37:05 


12 July 1995 


SHARP 


K 


c 


Ml 


9 


Kn Anon 2 


19:01:06.9 


-36:58:07 


12 July 1995 


SHARP 


K 


n 


GO (b) 


10 


V709 CrA 


19:01:34.9 


-37:00:57 


12 July 1995 


SHARP 


K 


w 


Kl 


11 


[MR81]HQr2 


19:01:41.6 


-36:59:53 


12 July 1995 


SHARP 


K 


c 


K8 


12 


V702 CrA 


19:02:02.0 


-37:07:44 


12 July 1995 


SHARP 


K 


w 


G5 


13 


CrAPMS 3 


19:02:22.1 


-36:55:41 


12 July 1995, 4 July 2001 


SHARP 


JHK 


w 


K2 


14 


[MR81] Ho- 14 


19:02:27.2 


-36:58:10 


12 July 1995, 2 July 2001 


SHARP 


JHK 


w 


MO-3 


15 


[GP75]RCrAe2 


19:01:27.2 


-36:59:09 


2 July 2001 


SHARP 


K 


T 


M3-5 


16 


[GP75]RCrAf2 


19:01:09.7 


-36:47:53 


2 July 2001 


SHARP 


K 


w 


K4 


17 


[GP75]RCrAn 


19:01:47.9 


-36:59:30 


2 July 2001 


SHARP 


K 


c 




18 


[MR81]Hq-12 


19:00:01.7 


-36:27:58 


6 July 2001 


SHARP 


K 


T 


M3-5 


19 


[MR81]Hq'13 


19:02:00.1 


-37:02:22 


6 July 2001 


SHARP 


K 


T 


M3-5 


20 


[MR81]Hq'15 


19:04:17.3 


-36:59:03 


6 July 2001 


SHARP 


K 


T 


M3-5 


21 


V721 CrA 


19:09:45.9 


-37:04:26 


2 July 2001 


SHARP 


K 


w 


K 


22 


[MR81]HQr 17 


19:10:43.4 


-36:59:09 


6 July 2001 


SHARP 


K 


T 


M3-5 


23 


CrAPMS 4 NW 


18:57:17.8 


-36:42:36 


4 June 2000 


ADONIS 


K 


w 


M0.5 


24 


CrAPMS 4 SE 


18:57:20.7 


-36:43:00 


4 June 2000 


ADONIS 


K 


w 


G5 


25 


CrAPMS 5 


18:58:01.7 


-36:53:45 


4 June 2000 


ADONIS 


K 


w 


K5 


26 


CrAPMS 6 


18:59:14.7 


-37:11:30 


5 July 2001 


SHARP 


K 


w 


M3-4 


27 


CrAPMS 8 


19:00:28.9 


-36:56:02 


2 July 2001 


SHARP 


K 


w 


M3 


28 


CrAPMS 9 


19:00:39.1 


-36:48:11 


5 July 2001 


SHARP 


K 


w 


M2 


29 


RXJ1855. 1-3754 


18:55:12.0 


-37:53:52 


5 June 2000 


ADONIS 


K 


w 


K3 


30 


RXJ1836.6-3451 


18:36:39.5 


-34:51:26 


6 June 2000 


ADONIS 


K 


w 


MO 


31 


RXJ1839.0-3726 


18:39:05.3 


-37:26:22 


5 June 2000 


ADONIS 


K 


w 


Kl 


32 


RXJ1840.8-3547 


18:40:53.3 


-35:46:45 


3 July 2001 


SHARP 


K 


w 


M4 


33 


RXJ1841.8-3525 


18:41:48.6 


-35:25:44 


5 June 2000 


ADONIS 


K 


w 


G7 


34 


RXJ1842.9-3532 


18:42:58.0 


-35:32:43 


6 June 2000 


ADONIS 


K 


c 


K2 


35 


RXJ1844.3-3541 


18:44:21.9 


-35:41:44 


5 June 2000 


ADONIS 


K 


w 


K5 


36 


RXJ1844.5-3723 


18:44:31.1 


-37:23:34 


6 June 2000 


ADONIS 


K 


w 


MO 


37 


RXJ1845.5-3750 


18:45:34.8 


-37:50:20 


6 July 2001 
5 June 2000 


SHARP 
ADONIS 


JHK 

K 


w 


G8 


38 


RXJ1846.7-3636 


18:46:45.6 


-36:36:18 


4 July 2001 

5 June 2000 


SHARP 
ADONIS 


JHK 

K 


w/w 


K6-7 (c) 


39 


RXJ1852.3-3700 


18:52:17.3 


-37:00:12 


6 June 2000 


ADONIS 


K 


c 


K3 


40 


RXJ1853. 1-3609 


18:53:06.0 


-36:10:23 


5 June 2000 


ADONIS 


K 


w 


K2 


41 


RXJ1856.6-3545 


18:56:44.0 


-35:45:32 


6 July 2001 
6 June 2000 


SHARP 
ADONIS 


JHK 
K 


w 


M2 


42 


RXJ1857.5-3732 


18:57:34.1 


-37:32:32 


2 July 2001 


SHARP 


K 


w/w 


M5-6 


43 


RXJ1901.4-3422 


19:01:28.7 


-34:22:36 


5 June 2000 


ADONIS 


K 


w 


F7 (d) 


44 


RXJ1901.6-3644 


19:01:40.5 


-36:44:32 


6 July 2001 
6 June 2000 


SHARP 
ADONIS 


JHK 

K 


w 


MO 


45 


RXJ1917.4-3756 


19:17:23.8 


-37:56:50 


3 July 2001 
5 June 2000 


SHARP 
ADONIS 


K 
K 


w 


K2 


46 


HR7169 


19:01:03.3 


-37:03:39 


1 July 2001 


SHARP 


K 


H 


B9 


47 


HR7170 


19:01:04.3 


-37:03:42 


1 July 2001 


SHARP 


K 


H 


B8 


48 


HD 176386 


19:01:38.9 


-36:53:27 


1 July 2001 


SHARP 


K 


H 


B9 


49 


SAO 210888 


19:04:44.4 


-36:50:41 


1 July 2001 


SHARP 


K 


H 


B9.5 



Notes: (a) c for cTTS, w for wTTS, T for TTS (unknown whether cTTS or wTTS), H for Herbig AeBe star, n for n on-TTS. (b) Non-mem bers 
JNeuhauser et al .1120001) . (c) Two separate TTS. (d) Foreground star at 65 ± 5 pc, not a member of the CrA association JNeuhauser et alj200Q) . (e) 
All other stars are members. 



In this program, the modulus of the complex visibility (i.e. the 
Fourier transform of the object brightness distribution) is deter- 
mined from power-spectrum analysis, the phase is computed us- 
ing the Knox-Thompson algorithm (iKnox & Thompson 1974<') 
and from the bispectrum dLohmann et al.lfl983h . Figure [U shows 



examples of reconstr ucted complex visib ilities. For a more de- 
tailed description see lKohler et al.l (|2000|) . 

If the object appeared unresolved, we computed the maxi- 
mum brightness ratio of a companion that could be hidden in 
the noise of the data. The principle was to determine how far 



[MR81]Hal4 
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[MRSllHo-lS 




CrA PMS 8 



RXJ1845.5-3750 







RXJ 1846.7-3636 NE 



i i^^fafei^ ■ 



RXJ1901.6-3644 







RXJ1917.4-3756 



HR7170 





Fig. 1. Complex visibilities of new binaries reconstructed from our SHARP I speckle-interferometric data. For each binary, the 
modulus of the complex visibility is shown on the left, and the phase computed using the Knox-Thompson algorithm on the right. 



the data deviated from the nominal result for a point source 
(modulus = 1, phase - 0), and then to compute the maximum 
brightness ratio of a companion that would be compatible with 
this amount of deviation. This was repeated for different posi- 
tion angles, and the maximum was used as upper limit for the 
brigh tness ratio of an undetected companion. See ' Leinert et al.l 



Prigl 
(1199 



1997ai) for a more detailed description of this procedure. 



If the complex visibility showed the fringe pattern typical of 
a binary, we computed a multidimensional least-squares fit us- 
ing the amoeba algorithm jPress et al.lll992l) . Our program tried 
to minimize the difference between modulus and phase com- 
puted from a model binary and the observational data by varying 
the separation, position angle, and brightness ratio of the model. 
This was necessary because the reconstructed images are a com- 
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CrA PMS 4NW 




Fig. 2. Like Fig. [T] but for new binaries discovered in our 
ADONIS data. 



plex function of the 2-dimensional separation vector and flux 
ratio that cannot be solved to compute the binary parameters di- 
rectly from the data. Fits to different subsets of the data yielded 
an estimate for the standard deviation of the binary parameters. 
We then subtracted the contribution of the companion from the 
images and applied the procedure described in the previous para- 
graph to find limits for the brightness of an undetected compan- 
ion. 

To derive the exact pixel scale and orientation of the detector, 
we took images of fields in the Orion Trapezium and the Galactic 
center during the observing campaigns at the NTT. The instru- 
mental positions of the stars i n the Trapezium were compared 
with the coordinates given in 'McCaughr ean & Stauffeij d 19941) 
by the astrometric software ASTR0M3- The pixel scale was 
49.1 ± 0.1 mas/pixel in 1995 and 49.5 + 0.2 mas /pixel in 2001. 
In 1995, the detector was rotated (90.5 + 0.1)° clockwise with 
respect to north on the sky, while it was rotated (90.2 + 0.2)° 
counterclockwise in 2001. We crosschecked the results with an 
astrometric calibration derived independently from images of the 



Fig. 3. Our K-band SHARP I images of wide binaries in Corona 
AustraHs. RXJ 1846.7-3636 NE-SW is not regarded as binary in 
this paper, because of the 6 arc sec upper separation limit chosen; 
RXJ 1846.7-3636 NE itself is a close binary shown in Fig.[I] 



see |http://www.starIink.rLac.uk/star/docs/sun5.htx/sun5.htniI| 



Galactic center (calibrated with iMenten et alJll997h and found 
that they agree within the uncertainties. 

During the ADONIS observations in June 2000, Orion 
could not be observed during the night. Instead, we used im- 
ages of fields near the Galactic center and positions given by 
Menten et al. (1997) for the astrometric calibration. The pixel 
scale was 49.6 ± 0.2 mas/pixel, and the detector was rotated 
(179.6 + 0.1)° clockwise. 
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Table 3. Companion candidates detected in this work 



Designation 


Date of 


Filter 


Separation 


Position 


Brightness Remarks 




Observation 




["] 


Angle [°] 


Ratio 


SCrA 


12 July 95 


K 


1.346 + 0.003 


149.3 + 0.1 


0.432 + 0.014 known 




6 July 01 


K 


1.330 + 0.004 


150.5 + 0.2 


0.401 +0.016 




6 July 01 


H 


1.340 + 0.009 


150.7 + 0.3 


0.356 + 0.012 




6 July 01 


J 


1.334 + 0.003 


150.3 + 0.2 


0.380 + 0.011 


VVCrA 


12 July 95 


K 


2.082 + 0.003 


43.5+0.1 


0.515+0.004 known 




6 July 01 


K 


2.079 + 0.007 


43.5 + 0.2 


0.358 + 0.013 




6 July 01 


H 


2.083 + 0.006 


43.7 + 0.2 


0.111 +0.008 




6 July 01 


J 


2.075+0.015 


43.8 + 0.2 


0.017+0.001 


CrAPMS 3 


12 July 95 


K 


4.403+0.018 


57.8 + 0.2 


0.146 + 0.003 known 




4 July 01 


K 


4.442 + 0.003 


57.5 + 0.2 


0.151 +0.005 




4 July 01 


H 


4.452 + 0.007 


57.5+0.2 


0.125+0.005 




4 July 01 


J 


4.455 + 0.009 


57.4 + 0.2 


0.106 + 0.005 


[MR81]Hq'14 


12 July 95 


K 


0.229 + 0.003 


341.0 + 0.4 


0.877 + 0.096 




2 July 01 


K 


0.216 + 0.003 


158.9 + 0.3 


0.891 +0.041 




2 July 01 


H 


0.214 + 0.005 


339.4 + 0.6 


0.860 + 0.062 




2 July 01 


J 


0.217 + 0.003 


338.7 + 0.6 


0.830 + 0.034 


[MR81] Ho- 15 


6 July 01 


K 


0.220 + 0.022 


278.1+2.8 


0.57 +0.285 


[MR81] Ha 17 


6 July 01 


K 


1.289 + 0.030 


46.4+1.9 


0.046 + 0.004 


CrAPMS 4 NW 


4 June 00 


K 


0.361+0.006 


261.7 + 0.2 


0.178 + 0.018 


CrAPMS 6 


5 July 01 


K 


2.363 + 0.009 


45.4 + 0.2 


0.835 + 0.043 known 


CrAPMS 8 


2 July 01 


K 


0.132 + 0.009 


110.4 + 4.3 


0.70 +0.300 


RXJ1844.3-3541 


5 June 00 


K 


0.227 + 0.003 


184.6 + 0.3 


0.396 + 0.023 


RXJ1845.5-3750 


5 June 00 


K 


0.904 + 0.003 


38.8 + 0.2 


0.602 + 0.007 




6 July 01 


K 


0.886 + 0.006 


38.6 + 0.2 


0.542 + 0.039 




6 July 01 


H 


0.898 + 0.004 


38.5 + 0.2 


0.568 + 0.027 




6 July 01 


J 


0.891+0.003 


38.5 + 0.4 


0.508 + 0.060 


RXJ 1846.7-3636 NE-SW 


5. June 00 


K 


7.847 + 0.003 


226.2 + 0.2 


0.371 +0.002 


NE-SW 


4 July 01 


K 


7.840 + 0.011 


226.0 + 0.2 


0.442 + 0.021 


NE-SW 


4 July 01 


J 


7.863 + 0.050 


226.0+1.0 


0.454 + 0.003 


RXJ1846.7-3636NEAB 


5 June 00 


K 


0.171 +0.004 


96.6 + 0.3 


0.577 + 0.031 


NEAB 


4 July 01 


K 


0.171+0.003 


99.9 + 0.3 


0.680 + 0.013 


NEAB 


4 July 01 


J 


0.168 + 0.003 


100.1+0.4 


0.620 + 0.015 


RXJ1853. 1-3609 


5 June 00 


K 


0.516 + 0.006 


88.7 + 0.6 


0.518 + 0.036 


RXJ1856.6-3545 


6 June 00 


K 


3.678 + 0.013 


88.2 + 0.2 


0.259 + 0.007 




6 July 01 


K 


3.644 + 0.003 


88.3 + 0.2 


0.226 + 0.003 




6 July 01 


H 


3.661+0.003 


88.3 + 0.2 


0.205+0.001 




6 July 01 


J 


3.670 + 0.006 


88.3 + 0.2 


0.223 + 0.003 


RXJ1857.5-3732 


2 July 01 


K 


2.861 + 0.003 


281.8 + 0.2 


0.516 + 0.011 known 


RXJ1901.6-3644 


6 June 00 


K 


0.942 + 0.003 


45.8 + 0.2 


0.458 + 0.010 




6 July 01 


K 


0.931+0.003 


45.3 + 0.4 


0.463 + 0.009 




6 July 01 


H 


0.935+0.014 


45.9 + 0.6 


0.469 + 0.020 




6 July 01 


J 


0.941+0.012 


46.2 + 0.5 


0.478+0.018 


RXJ1917.4-3756 


5 June 00 


K 


0.148 + 0.005 


198.7 + 1.4 


0.404 + 0.035 




3 July 01 


K 


0.138 + 0.003 


194.2 + 0.9 


0.350 + 0.024 


HR7170 


1 July 01 


K 


0.073 + 0.027 


98.6+13.4 


0.30 +0.153 


HD 176386 


1 July 01 


K 


4.104 + 0.004 


137.4 + 0.2 


0.128 + 0.013 known 



3. Results 

We have observed the sample of 49 stars listed in Table |2] Only 
a few known TTS members of CrA were not observed by us, 
namely the stars Patten Rlc, R17c, and R13a, as well as the M6- 
8 objects LS-CrA 1 & 2, Der iisl859, and CrA 432 & 444 , all 
discovered after our runs. See lNeuhauser & ForbrichI (l2008h for 
a full up-to-date list of all optically visible young members of 
R CrA. 

Of those 49 stars or systems, we now dismiss two ([MR 81] 
Ha 10 and Kn Anon 2) from the statistica l analysis, because 
they are nonmembers ( JNeuhauser et al.l2000l) . We also disregard 
RXJ 190 1.4-3422 from statistics, because it has a Hipparcos- 



measured distance of just 65 + 5 pc ; i.e. it is a foreground y oung 
star, but not a member of R CrA dNeuhauser et al.ll2005) . (We 
note that neither [MR 8 1 ] Ha 1 nor Kn Anon 2 nor RXJ 1 90 1 .4- 
3422 was found to be multiple.) Then, the stars RXJ1846.7-3636 
NW and SE are counted as two separate TTS, because their sep- 
aration is about 8 arcsec, more than the outer limit we adopt to 
avoid contamination by chance alignments (see Sect. 13. 2b . Those 
two objects were observed in one single observation, both within 
the SHARP-I field-of-view. The stars CrAPMS 4 NW and SE are 
also two separate TTS, because their separation is about one arc 
minute. We are left with 47 separately counted R CrA member 
systems observed, seven Herbig Ae/Be stars, and 40 TTS. In this 
sample of 47 member systems observed, we find 19 binary stars. 
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Results on separations, position angles (measured from N 
over E to S), and the corresponding flux ratios in the filters used 
are listed in Table [3] for all companions observed in this work, 
both previously known companions that we detected again and 
newly discovered companion candidates. The newly discovered 
binaries are shown in Figs.[T]|2l and[3] Figure|4]shows separation 
and brightness ratio of all companions. 
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Fig. 4. The results of our multiplicity survey in a plot of flux ra- 
tio or magnitude difference vs. binary star separation. The data 
points mark the detected companion stars, the thick line shows 
the average, and the thin line the worst sensitivity for undetected 
companions. The dashed vertical line at 0.13" shows the diffrac- 
tion limit for a 3.5 m telescope at K. This is the limit for unam- 
biguous identification of binary stars. The dashed horizontal line 
shows the completeness limit in flux ratio for the whole survey. 



3.1. Detection limits and completeness 

The modulus of the complex visibility of a binary is a cosine- 
shaped function (Fig.[TJ. If the separation of the binary is equal 
to the diffraction limit (0.13" for a 3.5 m telescope at K), ex- 
actly one period of the modulus of the visibility fits within the 
radius where the optical transfer function of the telescope is not 
zero. Under good circumstances, it is possible to discover bina- 
ries with even smaller separations, down to about half the diffrac- 
tion limitQ 

However, in these cases we can only detect the first mini- 
mum, but not the second maximum of the modulus of the visi- 
bility. Therefore, we cannot distinguish a close binary star with 
certainty from an elongated structure. Even more important, we 
cannot be sure that we find all companions at separations less 
than the diffraction limit. For these reasons, we limit ourselves to 
companions in the separation range between 0.13" and 6". The 
upper limit was chosen so that contamination with background 
stars has little effect (see Sect. 13.21 for a detailed discussion of 
this problem). 

Figure |4] shows not only the stars where we find compan- 
ions, but also the sensitivity of our survey, i. e. the maximum 
brightness ratio of a possible undetected companion as a func- 
tion of the separation. This sensitivity depends on factors like 
the atmospheric conditions at the time of the observations and 
the brightness of the target star We computed sensitivity limits 



^ We detected a close faint comp anion 0.073 arc sec east of HR 7170 
(late B-type SB according to Wils on & Jo^^l 19520 . a projected separa- 
tion of only 4 to 8 AU at the Hipparcos distance of 80 ± 20 pc. 



for each star in our survey, both stars where we find a compan- 
ion and stars where we do not find one. Fig.|4]shows the average 
sensitivity limit and the maximum of all limits, i.e. the envelope 
of the worst cases. 

Only the data of a few stars are too noisy to allow detec- 
tion of all companions down to magnitude differences to the pri- 
mary of 2.5 mag (corresponding to a brightness ratio of 1:10). 
These stars are [MR81] Ha 13, [MR81] Ha 15, CrAPMS 8, and 
RXJ1855. 1-3754. Based on the number of companions actually 
found, we expect less than 0. 1 additional companions above a 
brightness ratio of 1:10 at separations > 0.13". Therefore, we 
are confident we have found all companions with separations be- 
tween 0.13" and 6" and with a magnitude difference of less than 
2.5 mag. 

3.2. Chance alignments with background stars 

We observed most of our companion candidates only on one oc- 
casion; there is no way to tell from these data if two stars indeed 
form a gravitationally bound system or if they are simply two un- 
related stars that happen to be close to each other projected onto 
the plane of the sky. To confirm companionship, we would need 
a 2nd-epoch image some time later to show that both objects are 
co-moving; or even better, o ne would have to find curvature in 
the orbital motion (see e.g. Neuhauser et al.ll2008l) . To estimate 
the number of chance projections, we used the 2MASS point 
source catalog JSkrutskie et al.l2006l) . Within the area of our sur- 
vey (Right Ascension 279°. . . 294°, DecHnation -41° ... - 33°), 
the catalog contains 26180 sources brighter than K=12mag, the 
sensitivity limit of our imaging survey for close separations; 
2MASS is complete down to about K ^ 14 mag. This corre- 
sponds to a surface density of 218 sources per square degree, or 
1.7-10"^ sources per square arcsecond. Therefore, within 6" dis- 
tance from one of our targets, we expect 49-6^ -iTr- 1.7 -10"^ = 0.1 
field stars. We conclude that we can safely assume that all the 
companions we find are indeed gravitationally bound to their 
primary. 



4. Comparison with previous observations 

Among the 40 TTS observed, we detected the five known bi- 
naries again and discovered 12 new binaries. TY CrA BC was 
outside the field-of-view of SHARP, and TY CrA D is too faint 
to be detected reliably, although we see some indication of bi- 
narity in our speckle data. There are no known multiples among 
the three TTS not observed by us. There are no known triples or 
higher order multiples among these 40 TTS. 

We did not detect the com panion t o T Cr A fou nd in spectro- 
astrometric observations by iBailevI (Il998l) and iTakami et al.l 
62003i) . This co mpanion was also n ot det ected in the IR speckl e 
observations of iGhez et all ( Il997h and iLeinert et all ( Il997bh . 
The companion is only detected through its effect in Ha, but 
Bailey (1998) reports that the width and displacement of the Ha 
line is inconsistent with emission nebulosity or an outflow. They 
conclude that T CrA is a binary, and suggest that the companion 
might have strong Ha emission, but is too faint in the K-band to 
be detected by speckle interferometry. 

Among the 7 Herbig Ae/Be stars observed, we resolved HR 
7170 (new) and HD 176386 (known. lPommanget & Nvsll2002l) . 
Among all seven known Herbig Ae/Be stars in R CrA, there is 
one quadruple (TY CrA as eclipsing SB plus two resolved ob- 
jects), one triple (HR 7170 as SB plus one new resolved com- 
panion discovered in this work), and four more binaries (HR 
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7169 and R CrA as SBs, T CrA and HD 176386 as resolved 
binaries jj. Hence, among a total of 50 member systems (43 TTS 
plus 7 Herbig Ae/Be stars), there are one quadruple, one triple, 
and 21 binaries, i.e. 23 multiplefl Since the triple system con- 
tains two, and the quadruple system three companions, there are 
26 companion stars. In addition, among the five M6-8 type very 
low-mass stars and brown dwarfs known in R CrA, there are 
one to two binaries; however, we have left out the brown dwarfs 
and brown dwarf candidates in our statistics, because the brown 
dwarf completeness is very poor in CrA. 

Given 23 multiples among 50 systems, we have a multiplic- 
ity fraction (number of multiples divided by number of systems) 
of 46 + 10 %. Counting all binaries within triples and quadruples, 
26 companions in 50 systems give a companion-star frequency 
of 52 ± 10 %. If we restrict these numbers to those systems ob- 
served in this work and to those companions within the com- 
pleteness limits of this work (magnitude difference of lower than 
2.5 mag for separations between 0. 13 and 6 arc sec), we have 17 
binaries among the 47 systems observed here, i.e. a percentage 
of 36 ± 9 %. Since we do not find any triples or higher order 
multiples, the multiplicity and companion-star frequency are the 
same. 

Since we have observed young CrA members both on the 
dark R CrA cloud and (new ROSAT TTS) around the cloud, we 
can also compare the multiplicity on-cloud with off-cloud: We 
have nine binaries among 30 on-cloud members (a = 18:56 to 
19:24 and 6 = 38° to 36°), i.e. 30 + 10 %, and also eight binaries 
among 17 off^-cloud new ROSAT TTS, i.e. 47 + 17 %. While the 
binaries percentage off-cloud is higher than on-cloud, there is no 
statistically significant difference. 

We can now compare these numbers to other star-forming re- 
gions, which have been observed before with a similar technique 
and similar completness and sensitivity limits. 

4. 1 . Comparison to Main-Sequence stars 

The bina ry survey most comm o nly us ed for comparison is the 
work of iDuquennov & Mavoj (Il991l hereafter DM91), who 
studied a sample of 164 solar-type main-sequence stars in the 
solar neighborhood. They used spectroscopic observations, com- 
plemented by direct imaging; therefore they give the distribu- 
tion of periods of the binaries. Before we can compare this to 
our results, we have to convert it into a distribution of pro- 
j ected separations. We follow the method described in iKohlerl 
( 1200 lb . In short, we simulate 10 million artificial binaries with 
orbital elements distributed according to DM91. Then we com- 
pute the fraction of those binaries that could have been detected 
by our observations, i.e. those having projected separations be- 
tween 0'.'13 and 6" at the distance of CrA (130pc). The result is 
that 38.4 binaries out of the DM91 sample of 164 systems fall 
into this separation range, which corresponds to a companion- 
star frequency of (23.4 + 3.8)% (where the error is computed 
according to binomial statistics). 



^ The two SBs HR 7170 and HR 7169 form a common proper motion 
pair with only 13 arc sec separation, i.e. may form together a quintuple 
(including the newly resolved visual companion to HR 7170). 

■* Among the 13 newly discovered binaries, three are detected only by 
AO with ADONIS, stars not observed with SHARP I in speckle mode; 
five binaries are detected with both ADONIS and SHARP I speckle in- 
terferometry; five detected binaries were observed only with SHARP I 
speckle; among the new binaries found with SHARP I, five were de- 
tected only with speckle interferometry, not by simple shift+add imag- 
ing. 



Among the 47 systems in CrA, we find 17 binaries, which 
corresponds to 36.2+8.8 companions per 100 systems, 12.8+9.6 
more than among main-sequence stars. Thus, the companion- 
star frequency in our sample is higher by a factor of 1.5 + 0.4 
compared to main-sequence stars. Figure |5] shows the separation 
distribution of our sample and the sample of DM91. The over- 
abundance of binaries in CrA is somewhat more pronounced for 
larger separations, which might indicate that the peak of the sep- 
aration distribution is shifted to larger separations in CrA com- 
pared to main-sequence stars. However, due to the relatively 
small number of stars surveyed in CrA, this is not statistically 
significant. 

The TTS observed in our study have spectral types from 
G5 to M6 with luminosity classes III and IV and with ages on 
the order of one to a few Myrs, so th at they have masses be - 
tween ~ 0.1 and 1.5 M© according to iNeuhauser et al] ( |2000|) . 
These are somewhat lower masses than in the DM91 sample. 
Unfortunately, no multiplicity surve y of M-dwarfs comparable 
in size to DM91 has been published. JReid & GizisI (1 19971) stud- 
ied the largest sample so far, which contains 81 late K- or M- 
dwarfs. They find 22 companions in the separation range of 
our survey in CrA, which yields a companion-star frequency 
of (27 + 6) %. Within the errors, this is comparable to the 
companion-star frequency of DM91 in the s ame separation range 
(although the overall multiplicity found bv lReid & GizisI (1 19971) 
is lower than in DM91). We conclude that the companion-star 
frequency in CrA is also higher than among main-sequence M- 
dwarfs, by a factor of 1.3 + 0.3. 

If we assume that the separation distributions of our sample 
in CrA and the main-sequence sample of DM91 are the same, 
we can extrapolate the number of companions to all separations. 
Of DM9rs 101 companions, 38.4 fall on average into the sep- 
aration range surveyed by us. This means that the extrapolation 
factor is 101/38.4, which yields a total companion-star frequency 
of about (95 +23) %.This does «of imply that 95 % of the stars in 
CrA are multiple, but only that the average number of compan- 
ions per primary is about 0.95. Since, e.g., triple systems contain 
two companions and one primary, the fraction of multiples can 
be lower than 95 %. We do not know the ratio of binary to higher 
order multiples in CrA, therefore we have no way to extrapolate 
the total number of multiples. 

4.2. Comparison to otiner star-forming regions 

Table |4] lists the results of multiplicity surveys in several star- 
forming regions. Since the separation ranges where these sur- 
veys could detect companions are not always the same, it is 
not useful to compare the companion-star frequencies directly. 
Instead, it is standard practice to divide the companion-star fre- 
quency by the frequency of binaries among main-sequence stars 
(DM91) in the same separation range, and compare the resulting 
relative factors. 

In most T- and OB-associations, the companion-star fre- 
quency is found to be significantly higher than among main- 
sequence stars, by factors of 1.5 - 2 or even more. CrA shows 
the same high number of binaries and multiples, with a factor of 
1.5 + 0.4. 

If the separations of all the binaries found in two surveys are 
published, it is possible to count the companions in the separa- 
tion range common to both surveys and to compare the num- 
bers directly. We did this with our survey in CrA and the survey 
in Tau rus-Auriga by Leinert et al. ( 1993|) and[K6hler & Leinerj 
( Il998l) . The separation range in CrA (0'.'13 - 6") corresponds 
to 16.9 - 780 AU, or a.' 12 - 5'.'57 at the distance of Taurus- 
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Table 4. Multiplicity in different star-forming regions 



Region 



B :T: 

(a) 



Multiplicity 
[%] 



Companion-Star 
Frequency [%] 



relative 
to DM91 



Separation 
range ["] 



Ref. Remarks 



Tau-Aur 


60 


39 


3 


2 


42.3 ± 6.4 


49.0 ± 6.9 


1.9 ±0.3 


0.13^ 


-13 


(1) 


on-cloud 


Tau-Aur 


24 


21 








46.7 ± 10.2 


46.7 ± 10.2 


3.8 ±1.2 


0.09^ 


-2.5 


(2) 


on-cloud 


Tau-Aur 


40 


27 


2 


1 


42.9 ± 5.5 


48.6 ±8.3 


1.9 ±0.4 


0.13- 


-13 


(3) 


olF-cloud 


Lupus 


22 


10 


1 





33.3 ±10.1 


36.4 ±10.4 


1.4 ±0.4 


0.10- 


-12 


(4) 


on-cloud 


Cha 


19 


9 


1 





52.6 ± 16.6 


57.9 ± 17.4 


2.2 ±0.7 


0.10- 


-12 


(4) 


on-cloud 


Cha 


66 


11 








14.3 ±4.3 


14.3 ±4.3 


0.6 ±0.3 


0.13- 


-6 


(5) 


olF-cloud 


Oph-Sco 


13 


11 








45.8 ±13.8 


45.8 ±13.8 


2.9 ±0.9 


0.09- 


-2.5 


(2) 


on-cloud 


Oph 


114 


:42 


:2 


:0 


27.8 ± 4.2 


29.1 ±4.3 


1.2 ±0.3 


0.13- 


-6.4 


(6) 


on-cloud 


Sco-Cen 


59 


27 


2 





32.6 ±6.1 


35.2 ±6.3 


1.6 ±0.3 


0.13- 


-6 


(7) 


olF-cloud 



CrA 


30: 


17 : 0:0 


36.2 ± 8.8 


36.2 ±8.8 


1.5 ±0.4 


0.13-6 


(8) 


total 


CrA 


21 


:9:0:0 


30.0 ± 10.0 


30.0 ± 10.0 


1.3 ±0.5 


0.13-6 


(8) 


on-cloud 


CrA 


9 


8:0:0 


47.1 ±16.6 


47.1 ±16.6 


2.0 ±0.7 


0.13-6 


(8) 


off-cloud 



Remarks: (a) n umber of sing l es (S) to b inaries (B) to tripl es ( T) to quadruples (Q ). 

Refere nces : (DlLeinert et al.i(ll993h : (2) lGhez et alHl993h : (3) lK6hler & LeinertHl998h ; (4) lGhez et al1 ( ll997h : (5i lKohleil ( [200lh : (6) lRatzka et all 
( I2005h : (7) lK6hleretalJ l l200(]h :(8) this work. 
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Fig. 5. Companion-star frequency as a function of projected sep- 
aration. The histogram shows the result of this work, and the 
curve shows the distribution of separations for main-sequence 
binaries. 



Auriga . Th is range is fully covered b y the survey of' Leinert et al.l 
( Il993h and iKohler & Leinerj (Il998l) . They find 65 companions 
in this separation range, resulting in a companion-star frequency 
of 37.4 + 4.6 %. In CrA, we find 36.2 ± 8.8 %, which is almost 
identical. We conclude that the companion-star frequency in CrA 
is as high as in similar star-forming regions. 

5. Summary and conclusions 

We found 13 new binaries among young stars in the R CrA asso- 
ciation by simple speckle imaging, speckle interferometry, and 
AO imaging, all with 3.5-3.6 meter telescopes. 

The companion found next to [MR 81] Ha 17 should be 
mentioned explicitely (Fig. [3]l: It has a brightness ratio of only 
0.046 + 0.004 between the faint NE component and the bright 
SW component (Table |3]l. The unresolved system (i.e. t he bright 
component) has a spectral type of M3-5 ( lPattenlll998l) . Hence, 
the fainter object (K ^ 13.1 mag) may be a substellar compan- 
ion or a so-called infrared companion (e.g. with edge-on circum- 
stellar disk) or a background object. We continue to observe it in 



more detail (AO at larger telescopes and spectroscopy). All other 
companions found are almost certainly normal TTS. 

Among eight known cTTS, there are only two binaries 
(25 ±18 %), but among 28 known wTTS, there are 13 binaries 
found (46 + 13 %). No mixed systems were found, but for most 
companions and even some primaries, its not known whether 
they are wTTS or cTTS. The system [MR 81] Ho' 17 SW-NE 
may be mixed. 

The multiplicity in CrA is high and not significantly different 
from other similar star-forming regions like Sco-Cen or Tau-Aur 
The multiplicity among new ROSAT TTS around the dark cloud 
is high, just as it is in the case Tau-Aur There is no significant 
difference between on-cloud TTS and off-cloud TTS in CrA, but 
a tendency to more binaries among the off-cloud TTS. 

The multiplicity among the seven Herbig Ae/Be stars is high: 
one quadruple, one triple, and four binaries (from SBs to 6 arc- 
sec separation). If going out to 13 arc sec separation, we would 
even have one quintuple. However, the sample of known Herbig 
Ae/Be stars in CrA is too small to compare to other regions. All 
triples and quadruples known in CrA are hierarchical. 

For a comparison of the distribution of separations (or or- 
bital periods) found here in CrA with those of nearby solar-type, 
main-sequence stars (DM91), see Fig.|5] We can extrapolate to 
the total companion-star frequency (at any separation) in CrA 
being ~ 95 %, also similar to Sco-Cen and Tau-Aur. 
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